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Abstract 


This study examines the mechanisms of the Pacific decadal oscillation (PDO) in the 
GEOS-5 general circulation model. The model simulates a realistic PDO pattern that is 
resolved as the first empirical orthogonal function (EOF) of winter sea surface 
temperature (SST). The simulated PDO is primarily forced by Aleutian low through 
Ekman transport and surface fluxes, and shows a red spectrum without any preferred 
periodicity. This differs from the observations, which indicate a greater role of El Nino- 
Southem Oscillation (ENSO) forcing, and likely reflects the too short time scale of the 
simulated ENSO. The geostrophic transport in response to the Aleutian low is limited to 
the Kuroshio-Oyashio Extension, and is unlikely the main controlling factor in this 
model, although it reinforces the Ekman-induced SST anomalies. The delay between the 
Aleutian low and the PDO is relatively short (1 year) suggesting that the fast Ekman 
response (rather than Rossby wave propagation) sets the SST pattern immediately 
following an Aleutian low fluctuation. The atmospheric feedback (response to the SST) is 
only about 25% of the forcing and never evolves into an Aleutian low completely, instead 
projecting onto the North Pacific Oscillation (NPO), a meridional dipole in sea level 
pressure (SLP). The lack of preferred periodicity and weak atmospheric response both 
indicate a coupled oscillation is an unlikely mechanism for the PDO in this model. In 
agreement with recent studies, the NPO is correlated with the North Pacific Gyre 
Oscillation (NPGO), which is another leading EOF of the North Pacific SST. A possible 
connection between the PDO and the NPGO is discussed. 
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1. Introduction 


The dominant pattern of sea surface temperature (SST) variability in the extra-tropical 
Pacific is characterized by same-signed anomalies in the central and western parts of the 
basin and opposite signed anomalies along the west coast of the United States and the 
Gulf of Alaska, and is commonly referred as the Pacific decadal oscillation (PDO) or 
variability (Trenberth and Hurrell 1994; Mantua et al. 1997; Zhang et al. 1997). This 
pattern exhibits decadal to multi-decadal variability with marked ‘regime shifts’ around 
1925, 1947 and 1976 along with interannual variability that is largely in tune with the 
tropical El Nino-Southern Oscillation (ENSO) phenomenon. The decadal spectral peaks 
of the associated time series are broad and weak in the observational records, and 
therefore, the existence of a preferred oscillation is unclear (see, e.g., reviews by Miller 
and Schneider 2000; Liu et al. 2012). Nevertheless, pronounced low-frequency variability 
is apparent. A positive (negative) phase of the PDO is defined as the period when the 
eastern Pacific is anomalously warm (cool) and the central and the west Pacific is 
anomalously cool (warm) and is accompanied by negative (positive) sea level pressure 
(SLP) anomalies over the Aleutian Islands. 

The PDO appears to force atmospheric teleconnection patterns controlling climate 
variability in distant locations (Trenberth and Hurrell 1994; Mantua et al. 1997; Deser et 
al. 2004). For example, a positive phase of PDO coincides with enhanced wintertime 
precipitation over Alaska, southern United States and northern Mexico, reduced 
precipitation over much of the interior United States, and warm winter air temperature 
and reduced snowpack in the Pacific Northwest (Mantua et al. 1997). These 
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teleconnection patterns may be useful in seasonal and interannual prediction efforts in 
those regions. Additionally, understanding the relative role of natural decadal inodes like 
the PDO and the anthropogenic warming trend is important for near-tenn decadal 
prediction. 

A widely recognized model for the excitation of low frequency variability in extra- 
tropical oceans, in general, is that of the stochastic noise forcing (Hassehnann 1976), 
where SST variability is solely forced by the atmospheric ‘noise’. The ‘noise’ stands for 
atmospheric variability that is not forced by SST or other boundary states. In its simplest 
form, this model proposes that an atmospheric heat flux of white spectrum when coupled 
to a slab ocean model produces low frequency SST variability as the input heat flux is 
slowly damped due to the large heat capacity of the ocean. The resulting SST spectrum is 
red without any preferred periodicity. Extensions of this model include a propagative 
stochastic model where advection or wave propagation is taken into account (Frankignoul 
et al. 1997; Jin 1997; Saravanan and McWilliams 1997) and stochastically driven ocean 
dynamics model (Schneider et al. 2002). 

Another possibility is that the PDO is forced by ENSO through its atmospheric 
teleconnection to the North Pacific (Trenberth and Hurrell 1994; Zhang et al. 1997; 
Newman et al. 2003; Deser et al. 2004; Vimont 2005; Alexander and Scott 2008). This 
idea is supported by the observation that the two patterns are largely similar and the two 
modes underwent major decadal regime shifts in the observed records (Deser et al. 2004). 
The ENSO-forced atmospheric circulation patterns can generate SST anomalies in the 
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north-central Pacific on the seasonal timescale. These SST anomalies gain maximum 
amplitude in the spring following an ENSO mature phase, and can be stored beneath the 
mixed layer and resurfaced in the next winter via ‘reemergence’ (Alexander et al. 1999). 
The reemergence mechanism acts to enhance the persistence of SST anomalies from 
winter-to-winter thereby aiding in the reddening process (Deser et al. 2003). Thus the 
ENSO forced circulation anomalies can be considered as a ‘signal’ in the white noise- 
atmosphere that imparts temperature anomalies at quasi-regular intervals which has some 
persistence in the North Pacific Ocean mixed layer. Newman et al. (2003) showed that a 
first order autoregressive model (AR1) with an ENSO forcing forecasts the observed 
PDO with remarkable skill. Another viewpoint is that the PDO is due to the decadal 
modulation of ENSO, and therefore a statistical residue of ENSO on decadal timescales 
(Zhang et al. 1997; Vimont 2005). 

Latif and Barnett (1994) proposed that mid-latitude atmosphere-ocean interaction 
together with westward propagating oceanic Rossby waves could support a self-sustained 
decadal oscillation. Their proposal begins by assuming an SST anomaly in the Kuroshio- 
Oyashio extension (KOE) region of the subtropical gyre. This subsequently grows 
through Bjerknes-like ocean- atmosphere feedback and imparts a wind stress pattern in 
the central Pacific that in turn results in Ekman transport and Rossby wave propagation 
from the central Pacific to the western boundary. The resulting gyre anomalies eventually 
replace the original KOE SST anomalies with that of opposite sign. However, there is 
now a general consensus from observations and general circulation models (GCMs) that 
the central Pacific wind stress anomalies in fact lead the KOE SST and therefore are 
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unlikely due to an oceanic forcing (Schneider et al. 2002). This leads to a modified Latif- 
Barnett model in which a basin scale wind stress pattern with its maximum slightly 
shifted to the east of dateline forces oceanic Rossby waves in the central Pacific which 
generates KOE SST anomalies after 4-5 years (Deser et al. 1996; Miller et al. 1998; 
Schneider et al. 2002). The wind stress pattern is associated with the sea level pressure 
(SLP) anomalies of the Aleutian low. In this case, whether or not a closed oscillation 
occurs depends on if and how the atmosphere responds to the KOE SST. If the 
atmospheric response reverses the sign of the original wind stress curl anomalies, a 
coupled self-sustained oscillation may evolve. Results from the NCAR CCSM2.0 
(Known and Deser 2007) suggest that the PDO in that model is in agreement with the 
Aleutian low-induced Rossby wave mechanism. 

The length of observational records is a main constraint in evaluating these mechanisms, 
and it may be helpful to analyze long simulations of the atmosphere-ocean general 
circulation models (AOGCM), which are increasingly becoming computationally 
cheaper. In this study, we take advantage of a 350-year long simulation of the Goddard 
Earth Observing System (GEOS-5) AOGCM to examine the characteristics and 
mechanisms of the PDO. We will first examine the relative importance of the tropical 
versus extra-tropical forcing on the PDO - i.e., to what extent the ENSO-forced SST 
anomalies are reflected in the PDO. Secondly, we will examine the role of midlatitude 
atmospheric variability and the relative roles of surface heat flux, Ekman and geostrophic 
advection terms. In doing so, we hope to quantify the relative contributions to the PDO of 
a fast local response that occurs through surface heat flux and Ekman layer advection and 
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a delayed remote response at the KOE region, involving geostrophic adjustments which 
is indicative of the wind stress curl induced Rossby wave mechanism. 

While most previous studies of Pacific decadal variability have focused on the PDO, 
which is usually resolved as the first EOF of North Pacific SST variability, some recent 
studies have pointed out the importance of the second EOF (Bond et al. 2003), 
particularly in the years after the 1976 decadal shift. Di Lorenzo et al. (2008) found that 
the second EOF is connected to the eastern and central branches of the subtropical gyre 
and is often referred to as the North Pacific Gyre Oscillation (NPGO). In this study, we 
will briefly touch upon the NPGO in GEOS-5 and its potential connection with the PDO. 

Section 2 describes the model and observational data. Section 3 discusses the results, 
which include a brief summary of the mean climate (3.1), the spatial and temporal 
characteristics of the PDO and the NPGO (3.2), the tropical connection of the PDO (3.3), 
and extra-tropical atmospheric forcing and air-sea feedback (3.4). The summary and 
discussion is provided in Section 4. 

2. Model and data 

The atmospheric component of the GEOS-5 AOGCM was developed at NASA’s Global 
Modeling and Assimilation Office (Rienecker et al. 2008; Molod et al. 2012). Some key 
components of the model physics include the Relaxed Arakawa-Schubert scheme for 
convection (Moorthi and Suarez 1992), a Monin-Obukhov surface layer, which includes 
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the effects of a viscous sublayer for heat and moisture transport (Helfand and Schubert 
1995), and turbulence schemes by Lock (2000) and Louis and Geleyn (1982). The Ocean 
component is the Modular Ocean Model Version 4 (MOM4) developed by the 
Geophysical Fluid Dynamics Laboratory (Griffies et al. 2005), and the sea ice model is 
the Los Alamos sea ice model (CICE) (Hunke and Lipscomp 2008). The land surface 
model is a catchment-based scheme (Koster et al. 2000) that accounts for sub-grid scale 
heterogeneity in surface moisture. 

The resolution of the atmospheric model is 2.5° longitude by 2° latitude with 72 vertical 
layers and that of the ocean is 1° in longitude and latitude telescoping to 1/3° meridional 
spacing near the equator, with 50 vertical levels. The ocean vertical grid spacing is a 
constant 10 m over the top 225 m. The current simulation was started from an earlier run 
where the atmospheric model was initialized from an uncoupled run forced by observed 
SST. The ocean model was initialized with climatological temperature and salinity from 
Levitus and Boyer (1994a;b). Green house gas concentrations are fixed at 1950’s level. 
In this study, we examine the last 350 years from a 407-year long single coupled run. 
The first 57 years are omitted as model spin-up, as determined by visual examination of 
the globally averaged annual mean SST. 

The model results are compared with observed SST product of Rayner et al. (2003) for 
the period 1870-2011 and sea level pressure (SLP) for the period 1871-2010 from the 
NOAA 20 th century reanalysis (Compo et al. 2011). A global warming signal was 
removed from the observed SST by linear de-trending. Oceanic variables from the model 
are also de -trended to eliminate any model drift. Monthly or seasonal anomalies are 
obtained by removing their respective long-term mean values. A 6-year low-pass Lanczos 
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filter (Duchon, 1979) is utilized in certain analyses to better extract the decadal signals. 
Diagnostic tools include Empirical Orthogonal function (EOF) decomposition, spectral 
analysis, and lag correlation or regression. The statistical significance of the spectral 
peaks is estimated with a Chi-squared test with respect to a null hypothesis based on an 
AR1 red noise process. A two-sided t-test is used to assess significance of temporal 
correlation coefficients. The effective number of degrees of freedom for the t-test is 

calculated as, N' — N [LZhIiz j where N is the number of time steps in the data, N' is 

\l+T lx r ly J' F 

the effective number of degrees of freedom and r lx and r ly are the lag- 1 autocorrelation 

coefficients for time series x and y, respectively (Bretherton et al. 1999). 

3. Results 

3.1. Mean climate 

The GEOS-5 AOGCM produces a stable, realistic mean state, and realistic major modes 
of variability (Vikhliaev et al. 2011). This is confirmed in the current 350-year long 
simulation. For example, Fig. 1 shows that the model captures the climatological features 
of the tropical SST such as the Pacific warm pool and cold tongue fairly well. Some 
errors in the model are common to many current climate models and they include the 
westward extension of the cold tongue, which often results in a ‘double ITCZ’ pattern in 
precipitation, warm bias along the continental coasts, and a meridionally confined ENSO 
that extends too far to the west. The double ITCZ problem is often related to errors in 
atmosphere-ocean feedbacks in the tropical Pacific (Sun et al. 2003). The wann bias over 
the western side of the continents is likely due to the lack of low clouds in this model 
(Molod et al. 2012) or coastal upwelling biases prevalent in AOGCMs in general (e.g., 
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Griffies et al. 2009). The enhanced variance in the tropical Indian Ocean implies an 
unrealistically strong forcing from the Pacific, which may be due to the westward 
extension of the ENSO. 

3.2. Low frequency variability in North Pacific SST 

In this section, we describe the leading patterns of variance in the North Pacific SST that 
are isolated by an EOF analysis. These are the PDO and the NPGO, and we will refer to 
them as the leading SST modes throughout the remainder of the paper. The EOF analysis 
is applied to unfiltered December through February (DJF) seasonal mean anomalies of 
SST over the region 120°E-100°W, 20°N-60°N. The corresponding global spatial 
patterns are obtained by regressing the normalized principal components (PCs) on the 
DJF SST anomalies. Figure 2 shows these regressed fields based on 350 years of the 
model run and 141 years (1871-2011) of the Hadley Center SST. The first EOF captures 
the familiar horseshoe pattern of the PDO and explains 29% of the normalized variance 
in the model and 30% in the observations out of their total variance values which are 0.36 
and 0.21 (°C) 2 , respectively (Figs. 2a, c). 

Recently, there has been some interest in another leading mode of North Pacific SST 
variability, including its timescale, mechanisms and its role in the Pacific interannual- 
decadal variability (Bond et al. 2003; Di Forenzo et al. 2008). Found here as EOF2 in the 
model and EOF3 in the observations, it depicts a nearly north-south dipole structure, 
which is similar to the SST pattern in Bond et al. (2003) (Victoria mode) and the North 
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Pacific Gyre Oscillation (NPGO) in Di Lorenzo et al. (2008). We note that the ra nk of 
this mode appears sensitive to the observational period; while the NPGO is resolved as 
EOF2 for the period 1950-2011, its variance is reduced over the period 1871-2011. 

We focus on the EOFs for the DJF season in order to better capture the PDO-ENSO 
relation in the model. We note that although an EOF analysis of the annual mean SST 
yields similar patterns and variance decomposition, the simulated ENSO-PDO covariance 
is absent in that case: the regression coefficients in the ENSO region are less than 0.05. 
The regression values for DJF are, however, still small (~0.2) compared with those for 
the observations (-0.6) (Fig. 2c) and we will discuss this point further in the following 
sections. 

In Figure 2, there are two major differences between the model and observations. Firstly, 
the model overestimates the variance of both poles of the PDO mode. The negative 
anomalies in the western part of the basin (at about 40°N) and the positive anomalies 
over the Gulf of Alaska are nearly double that of their observational counterparts. The 
variance is particularly large in the western north Pacific region where the influence of 
the Kuroshio-Oyashio currents is important. The PDO pattern has potentially two centers 
of action, one in the western basin and another in the central part of the basin, the latter 
being related to ENSO (Deser and Blackmon 1993). In the observations, the central 
Pacific maximum appears dominant. In the model, this distinction is less clear and 
appears to be related to the nature of the ENSO, which has a pronounced biennial 
spectrum (figure not shown). The dominance of the KOE variance has also been seen in 
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different versions of the NCAR Community Climate System Model (CCSM), (Kwon and 
Deser, 2007; Deser et al. 2012) as well as in many models that participated in the 
Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment (AR4) (Furtado 
et al. 2011). 

Secondly, as mentioned before, the tropical signal is weak in the model. Compared to a 
broad ENSO-like pattern in the eastern Pacific found in the observations, the model 
shows weaker signals slightly shifted to the west. The apparent weak relation with ENSO 
is evident in the temporal correlation between PCI and Nino3: 0.17 in the seasonal DJF 
unsmoothed data and 0.34 when both time series are subjected to a 5-year running mean. 
(The Nino3 index is defined as the SST anomalies averaged over 150°W-90°W, 5°S- 
5°N). The corresponding value in the observations is 0.57 in the unfiltered data for the 
period 1870-2011. Reasons for this general weakness of the ENSO-PDO relationship will 
be examined in Section 3.3. 

The power spectrum of the PDO for the entire simulation (350 years) has three peaks 
around 31,18 and 7 years respectively, which are not statistically significant at the 10% 
level (figure not shown). Similarly, the NPGO has spectral peaks around 15 and 7 years 
that are not statistically significant (figure not shown). The lack of significance is even 
more evident when spectral power is computed independently for two or more segments 
of the data and averaged. What is shown in Fig. 3 is such an average over two 141-year 
segments, which shows a red spectrum with small increases in power near 28 and the 9-5 
year band. The results are similar when averaging is done over three 100-year segments. 
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Notice that no distinct peaks are present either in the PDO or in the NPGO, which is not 
inconsistent with the observations (Fig. 4) although questions remain about the accuracy 
of the observational records. 

This is different from the finding of Kwon and Deser (2007), where they found a 
statistically significant 16-year spectral peak in a 650-year long simulation of the 
CCSM2. A distinct spectral peak implies a preferred oscillation that could arise from 
either atmosphere-ocean feedback or from resonance. In CCSM2, the basin-wide wind 
stress curl forcing excites Rossby waves in the central Pacific and results in a delayed 
response in the surface temperature in the western boundary that in turn propagates 
eastward. These processes, wave propagation and advection combined, fonn an 8-year 
half cycle and closure of this cycle depends whether and how the atmosphere responds to 
the KOE SST. In CCSM2, the atmospheric response (in wind stress curl) to KOE SST is 
about 35% weaker than the original forcing; yet, the response is of opposite sign to the 
forcing, a necessary condition for establishing a coupled oscillation with a 16-year cycle. 
The fact that there is no distinct spectral peak in the GEOS-5 diminishes the prospects of 
a coupled self-sustained oscillation in this model. But that does not rule out the Rossby 
wave mechanism in this model - for example if the atmospheric response to the KOE 
SST anomalies is insufficient to generate a closed loop. Wind stress curl forced Rossby 
waves may still be present and influence the KOE SST anomalies and therefore the 
objectives of the study remain valid. 
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In the remaining sections of the paper, we mainly focus on the model data. The main 
reason for this is that the accuracy of the observational records before 1950 may be 
questionable and statistical calculations (such as regression) may not provide a 
benchmark to make meaningful comparisons against the model. Instead, the 
understanding gained from case studies such as the decadal shift of the 1976 (a phase 
change where data are more reliable) may be more useful and would help us to better 
evaluate the model. The wind-stress curl forced Rossby waves is a mechanism proposed 
to explain the 1976 phase shift (Deser et al. 1999; Seager et al. 2001) and we will take 
this as an observational benchmark in this study. 

3.3. Tropical connection 

As mentioned before, the observed relatively strong relation between the PDO and ENSO 
is not captured in the model. The simultaneous correlation between the PCI of SST and 
the Nino3 index is 0.17 in the unsmoothed DJF seasonal data in contrast to 0.45 in the 
observations. We speculate that the ENSO-induced SST fluctuations do not project onto 
the PDO in the simulation. A mechanism by which ENSO can generate low frequency 
fluctuations involves persistence of temperature anomalies in the mixed layer for a few 
years aided by ‘reemergence’. This may operate in nature where the ocean has a chance 
to retain an initial anomaly for 4-5 years between two consecutive ENSO events. In the 
model, however, the time scale of ENSO is too short (dominant period of 2-3 years) such 
that a positive ENSO event in one winter is often immediately followed by a negative 
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event in the next winter and therefore is unlikely to force low frequency variability in the 
North Pacific. 

The ENSO-forced signals in the North Pacific can also be examined in terms of Nino3- 
based monthly lagged composites of SST and SLP (Fig. 5). A wann phase is defined 
when DJF Nino3 exceeds +1 standard deviation and vice versa. The circulation response 
consisting of negative SLP anomalies during the warm phase and positive SLP anomalies 
during the cold phase is in agreement with observations. Several studies have examined 
the effect of this atmospheric teleconnection on the North Pacific SST (see, Alexander et 
al. 2002 for a review). During a warm ENSO phase, a negative SLP anomaly in the North 
Pacific results in anomalous northwesterlies and cold dry air advection to the 
northwestern Pacific, anomalous southerlies and warm moist air advection to the West 
coast of the North America and anomalous westerlies and increased total westerly winds 
in the central Pacific. The resulting surface heat flux and Ekman advection anomalies 
cool the ocean west of 150°W and warm the ocean on the eastern side. During a cold 
ENSO phase, these anomalies are reversed. Unlike the response during warm ENSO 
events, the response in SST for cold ENSO events is weak and incoherent in the model. 
In the cold composites, negative SST values occur in the North Pacific as late as 
November, suggesting that negative anomalies from a previous warm ENSO persist into 
the following fall. This provides an unfavorable condition for circulation anomalies 
associated with a cold ENSO phase to generate warm SST anomalies. This problem is 
clearer when composites are extended into the preceding fall and the following summer 
(figure not shown). During summer, the tropical Pacific reverses sign abruptly while the 
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extra-tropical SST anomalies retain the sign gained in the previous winter. This results in 
an erroneous Pacific SST pattern consisting of a cold tropics and subtropics or a warm 
tropics and subtropics. Interestingly, this behavior is obvious only in the cold composites. 
Studies that examined the asymmetric response to ENSO show that, it is the cold events 
that force larger circulation and SST anomalies than wann events (Hoerling et al. 1997; 
An et al. 2007), while our model composites show the opposite. The reasons for this 
asymmetry are not immediately clear from the present analysis. Nevertheless, the 
composite analysis verifies our earlier conclusion in Section 3.2 that the tropically forced 
SST anomalies in the model occur in the central Pacific slightly south of 40°N, which is 
in agreement with observations, while the PDO variance (see, Fig. 2) is maximum over 
the western basin. We summarize therefore that the influence of ENSO may in fact be 
contributing to interannual fluctuations in North Pacific SST, but not necessarily 
projected onto the decadal timescales studied here. 

3.4. Extra-tropical atmospheric forcing and feedback 

We now move on to looking more directly at the role of extra-tropical atmospheric 
forcing of the PDO. Several studies have pointed out the correlation between Aleutian 
low variability (often depicted by the North Pacific Index or NPI) and the PDO (Deser et 
al. 1999; Kwon and Deser, 2007; Schneider and Comuelle, 2005)). The NPI is an area- 
average of SLP anomalies over the region 30°N-65°N, 160°E-140°W as defined in 
Trenberth and Hurrell (1994), and captures a monopole pattern similar to the first EOF of 
winter SLP anomalies. However, we know that the influence of ENSO on the Aleutian 
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low variability is not negligible as the two indices (NPI and Nino3) are correlated at 0.5, 
indicating that the ENSO accounts for about 25% of the Aleutian low variance. 
Additionally, the ENSO composites of SLP (Fig. 5) indicate an atmospheric response 
over the Aleutian low region. Therefore, the Aleutian low variability is likely a 
combination of internal atmospheric and the ENSO-forced variability and in the 
following analysis we attempt to separate these two components by linear regression. 

a) AR1 forecast: 

We will first examine the influence of NPI by adopting a simple AR1 model with an 
additional forcing term to represent the NPI, following Newman et al. (2003). 
Accordingly, the PDO is represented as follows: 

P(t) = «P(t-l)+ 0S(t)+ 77(t), (1) 

where P(t) is the PDO forecast at time t, P(t — 1) is the PDO one time step earlier, S is 
an index to represent external forcing (e.g., NPI) and r] is an uncorrelated noise tenn. The 
coefficients fJ and oc are obtained, respectively, by regressing the PDO on the NPI time 
series and then auto-regressing the residual PDO times series (that is obtained after 
removing the NPI) with a lag of one year. Forecasts are made for one time step (or 1 
year) at a time and then reinitialized and run for the next time step and so on. To separate 
the ENSO-forced part of the NPI, a revised NPI (NPI*) is defined, which is obtained by 
removing the ENSO signal by linear regression with the Nino3. The NPI* is an 
approximation to the inherent (internal) extra-tropical atmospheric variability. We make 
two forecasts, first with the total NPI, and then with the NPI*. 
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Figure 6 shows scatter plots of the PDO versus the NPI and the NPI* (Figs 6a, b) and the 
AR1 forecasts of the PDO when run with NPI and NPI* (Figs. 6c, d). The correlation 
between the PDO and the NPI (i.e., /? in (1)) remains unchanged whether or not ENSO is 
present, and therefore, the skill of the forecast. This suggests that the ENSO-forced part 
of the NPI is not the key forcing of the PDO and that the ENSO-forcing is likely limited 
to the high frequency variability (1-2 years) of the NPI. These results are again consistent 
with our previous finding that the role of the ENSO on the PDO is negligible in this 
model. This finding is however different from those based on observations. For example, 
Schneider and Comuelle (2005) found that by removing ENSO from observed North 
Pacific SLP anomalies, the skill of the AR1 reconstruction of the PDO was reduced from 
~0.8 to -0.55. They identified a region, where intrinsic midlatitude SLP variability is 
important for the PDO reconstruction, which is approximately collocated with the NPI 
region of Trenberth and Hurrell (1994), used in this study. 

b) Atmospheric-Ocean feedback 

In view of the simultaneous covariance between the NPI/Aleutian low and the PDO, we 
next examine lead/lag relationships between the atmosphere and SST in order to 
understand forcing and response between atmosphere and SST. We remind the reader 
that the original Latif-Barnett hypothesis (Latif and Barnett, 1994) called for an oceanic 
origin of the SST in which case one would expect NPI anomalies to lag the PDO. There 
is however a general agreement in recent studies that the NPI leads the PDO (e.g., Kwon 
and Deser, 2007; Deser et al. 2012). Another recent finding is that the second mode of the 
SST (NPGO) is correlated with the North Pacific Oscillation (NPO) (Chhak et al. 2009). 
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The NPO is an SLP pattern identified by Walker and Bliss (1932) that has a meridional 
dipole structure, with opposite signed anomalies in Alaska and Hawaii. 

We first isolate the Aleutian low and NPO patterns by computing EOFs of DJF SLP and 
wind stress curl over the region 120°E-100°W, 20°N-60°N (Fig. 7). Here, instead of 
using the NPI index, we use the EOF analysis to identify the Aleutian low since this 
enables a cleaner separation of the Aleutian low and the NPO. The EOFs of the SLP and 
the wind stress curl essentially capture the same modes: the temporal correlation between 
the corresponding PCs are greater than 0.9. The first mode (Figs. 7a, c) shows the 
Aleutian low and its signature in wind stress curl, the temporal variability of which is 
essentially the NPI: PCI of the wind stress curl is correlated with the NPI with a value of 
0.97. The second mode captures the characteristic meridional dipole structure of the 
NPO. Consistent with internally forced midlatitude atmospheric behavior, (e.g., Wallace 
and Gutzler 1981; Esbensen 1984) the temporal variability of these modes at interannual 
time scales shows a white spectrum without any preferred periodicity (figure not shown). 
In the following analyses, we will use PCs 1 and 2 of the wind stress curl to represent the 
Aleutian low and NPO, respectively. 

Figure 8 shows lead/lag correlations between the PCs of the wind stress curl and the SST. 
All time series are subjected to a 6-year low-pass filter to better capture decadal-scale 
variability. First of all, considering the relation between the Aleutian low and the PDO, 
the correlation is at maximum (~0.7) when the atmospheric mode leads the SST mode by 
1 year, suggesting an atmospheric forcing of the SST (Fig. 8, bottom panel). Similarly, 


18 


460 

461 

462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 


the NPO appears to force the NPGO with a maximum correlation of ~0.4 at lag -1. The 1- 
year lag between NPO and the NPGO is in agreement with previous studies (Chhak et al. 
2009). In the case of the PDO, the 1-year lag is too short to invoke the wind stress curl 
induced Rossby wave mechanism, in which Rossby wave propagation from the central 
Pacific to the western boundary sets the delay time between the Aleutian low and the 
PDO SST anomaly. The short time lag obtained here is not in agreement with the results 
based on CCSM (Kwon and Deser, 2007; Deser et al. 2012), where lags of 2-4 years 
were found that are consistent with a Rossby wave mechanism. Rossby wave phase 
speed estimates for 40°N suggest 3-5 years are needed to propagate from 160°E to the 
western boundary (Chelton and Schlax, 1996). 

Secondly, we notice that there are weak correlations between the PDO and the NPGO 
when the PDO leads the NPGO by 3 years and similarly between the Aleutian low and 
the NPO when the former leads by 3 years (Fig. 8, top panel). These values are about 0.4 
and 0.2, respectively and they are marginally significant at the 5% level. A similar weak 
relationship (statistically significant at 5% significance level) can be found between the 
PDO and the NPO when the fonner leads the latter by 2 years, suggesting that the 
atmospheric response of the PDO projects onto the NPO. In view of the above 
correlations, we propose a connection between the PDO and the NPGO, which involves 
the following processes in the order they are listed: 1) an initial Aleutian low anomaly 
forces the PDO, 2) the atmosphere responds to the PDO which results in an NPO-like 
pattern, and 3) the NPO forces the NPGO SST pattern. 
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To further address the atmospheric forcing and response, lead lag regression and 
correlations are computed for SLP and wind stress curl data with respect to the PDO and 
the NPGO (Figs. 9, 10). Regression with SST is shown (left column in Figs 9, 10) to 
understand the progression of the EOF patterns. In the regressions with atmospheric 
variables, a positive lag means the atmosphere lags the ocean and therefore the regressed 
fields represent an atmospheric response and similarly a negative lag suggests 
atmospheric forcing. Only lags from -2 to +3 years are examined, as they appear relevant 
based on Fig. 8. In the case of PDO (Fig. 9), lags -2 to 0 show SLP and wind stress curl 
patterns that represent the Aleutian low, similar to the EOFs 1 in Fig. 7. The atmospheric 
response to the PDO (lags +1 to 3) is of opposite sign to the forcing, and forms over the 
western boundary, expanding eastward, south of 40°N. Although the response is of 
opposite sign to the forcing, it does not develop into the monopole structure of the 
Aleutian low, rather it bears similarities to the NPO pattern (see, Figs, 7c, d). Figure 10 
shows that an NPO-like structure in SLP and wind stress curl leads the NPGO, as one 
would expect from the correlation curve in Fig. 8. Also, notice that at lag -2 the SST and 
SLP patterns are similar to the PDO where the SST has maximum amplitude at the 
western boundary at 40°N, and the SLP has the Aleutian low structure. This is again 
consistent with Fig. 8, where we showed that the PDO leads the NPGO. 

In summary, the lead/lag regression analysis conveys the following. 1) It is unlikely that 
the atmospheric response to the PDO matures into the original forcing pattern of the 
Aleutian low. Instead, the response appears to project onto a meridional dipole in SLP, 
which is similar to the NPO pattern. 2) And therefore, the atmospheric response to the 
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PDO may help force the NPGO. 3) The atmospheric response is generally weak - a rough 
estimate from the regression analysis is that the response is 25% of the forcing. 


c) Oceanic response to the Aleutian low 

A key question now is how the SST anomalies are generated, i.e., whether they are due to 
surface heat fluxes and Ekman layer advection, processes that are fast and local, or by 
wind-stress curl induced Rossby waves (e.g., Deser et al. 1999; Seager et al. 2001; Kwon 
and Deser, 2007), a delayed response to the Aleutian low in the KOE region, which can 
take 3-5 years to generate SST anomalies in the KOE region. To separate these 
processes, we consider the following familiar approximation for the mixed layer heat 
budget (see, review by Frankignoul 1985) 


where T is the mixed layer temperature or equivalently SST, Q is the net surface heat 
flux, p 0 is the density of sea water (1025 kgm' 3 ), C p is the specific heat capacity of sea 


Ekman and geostrophic velocity, respectively. The rest of the terms form the residual, R 
that include diffusion and entrainment. The Ekman transport is computed from wind 
stress as, V E = kxr/ p 0 fh, and the surface geostrophic velocity is obtained from sea 
surface height as, V g — —kxVg(g / /), where r is the surface horizontal wind stress 
vector, Tj is the sea surface height and g and / are gravity and Coriolis parameter, 
respectively. Total fields of surface currents are obtained from the DJF mean wind stress 
and sea surface height from which DJF anomalies are calculated by subtracting the 
climatology. The mixed layer depth in the Ekman current calculation is the climatological 


dT 

dt 



( 2 ) 


water (3986 Ihg^K' 1 ), h is the mixed layer depth and V E and V g are horizontal vectors of 
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mean 2D field, which ranges from 80- 140m west of dateline between 30°N-50°N and is 
less than or equal to 60m over the rest of the basin. 

We will examine regressed fields of net heat flux, horizontal Ekman and geostrophic 
temperature advection (from (2)) with respect to an index of the PDO. It is not yet clear 
whether the NPGO is completely independent of the PDO and therefore, in this section 
we focus only on the PDO and its association with the KOE SST anomalies. For this 
purpose, we define an area-averaged SST index over the KOE region (140°E-180°E, 
35°N-45°N) (KOE index, hereafter), instead of using the PC time series of the PDO. 

Figure 1 1 shows regressed fields of the net surface heat flux with respect to the KOE 
index. The atmosphere leads or is simultaneous with the SST at negative and zero lags. 
At negative or zero lags, a positive regression coefficient implies either 1) heat flux is 
positive (into the ocean) and consequently SST increases or 2) heat flux is negative (into 
the atmosphere) and ocean surface cools in response. Positive values south of the KOE 
region and over the central parts of the basin, therefore suggests that those are the regions 
where the surface heat flux forces SST anomalies. On the other hand, negative values in 
the KOE region (40°N-50°N) implies that heat flux is out of the ocean as SST increases, 
suggesting that oceanic processes are important in controlling the SST, and surface heat 
flux works to only damp the SST. 

Figure 12 shows the regressed fields of anomalous Ekman and geostrophic advection of 
temperature with respect to the KOE index. Anomalous advection terms are computed 
from (2) by subtracting the climatological mean to obtain V ■ VT' , V' ■ VT, and V' ■ VT', 
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which are advection of anomalous temperature by mean currents, advection of mean 
temperature by anomalous currents and advection of anomalous temperature by 
anomalous currents, where V — (V E , F g ). The ‘overbar’ denotes climatology and ‘prime’ 
denotes deviations from climatology. The sum of the three terms forms the total 
anomalous advection. In Fig. 12, negative values denote areas where heat transport 
contributes to a positive KOE index. The Ekman temperature advection terms appear to 
determine the SST over most of the PDO pattern, especially east of 160°E. This is clear 
in panel (b) that shows the product of the anomalous Ekman current and the mean 
temperature gradient. The panel (a) represents advection of Ekman current across 
anomalous temperature gradients, which is dominant west of 160°E and between 38°N- 
42°N. Here, the anomalous temperature gradient is not necessarily only due to Ekman 
transport; it could be due to surface heat fluxes, geostrophic transport or other mixed 
layer processes that are not considered here. From Fig. 1 1, we see that surface heat flux 
term is important in the western boundary in areas south of 40°N, which partially 
overlaps our area of focus in panel (a). The anomalous geostrophic advection (panel h), 
on the other hand, is limited to a narrow strip above 42°N, west of 160°E, and it is 
unlikely a major contributor for panel (a). Therefore, the anomalous temperature 
gradients between 38°N-42°N that contribute to the term in panel (a), is unlikely due to 
geostrophic transport. The surface geostrophic advection represents the delayed 
response at the KOE, which can be manifested by Rossby waves from the east either 
resulting in the strengthening/weakening of the gyre (Latif and Barnett 1994) or by the 
meridional shift at the confluence of the subtropical and sub polar gyres (Seager et al. 
2001). Figure 12 indicates that in GEOS-5, the Ekman transport dominates over the slow 
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geostrophic advection. This is consistent with the short lag (1 year as opposed to 3-4 
years) between the Aleutian low and the PDO, shown in Fig. 8. 

In agreement with previous studies (Seager et al. 2001; Kwon and Deser 2007, among 
others), we find that the geostrophic anomalies are of the same sign as the Ekman 
anomalies. The phase reversal of the KOE SST, therefore, has its origins in the 
atmosphere, as proposed by recent studies as opposed to the ocean, which was the earlier 
proposal by Latif and Barnett (1994). The lead lag regressions in Fig. 9 indicate that the 
atmospheric response to the SST is rather weak and does not project on to the Aleutian 
low. This leaves us with the possibility that the PDO in this model is mostly an oceanic 
surface layer response to the Aleutian low and is unlikely a coupled ocean-atmosphere 
oscillation. 

4. Summary and discussion 

In this study, we examined the characteristics and mechanisms of Pacific decadal 
variability in a 350-year long simulation of the GEOS-5 atmosphere-ocean general 
circulation model. 

The GEOS-5 simulates a realistic pattern of decadal variability characterized by a PDO- 
like horseshoe pattern in the SST. The time series of this pattern does not have a 
statistically significant preferred oscillation, and can instead be best characterized by a 
red noise process. The PDO in this model is likely primarily forced by midlatitude 
atmospheric noise and is little influenced by tropical SSTs associated with the ENSO 
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phenomenon. The lack of a substantial role of ENSO may be related to the unrealistically 
short time scales of ENSO in the model of about 2-3 years, which is unlikely to be 
influential on the decadal timescale of the North Pacific SST. The lack of a substantial 
tropical connection must be a model bias since observational studies suggest ENSO as 
one of the triggering mechanisms for the PDO (Newman et al. 2003; Schneider and 
Cornuelle, 2005). 

The Aleutian low and the associated basin-wide wind stress curl pattern are correlated 
with the PDO at ~0.7 when the former leads the PDO by 1 year. The 1-year delay 
suggests that the local fast response of the ocean by surface heat fluxes and Ekman 
transport takes precedence over any delayed remote response at the KOE region. The 
E km an transport is dominant over most of the central and western parts of the PDO, 
while the influence of the geostrophic transport is limited to a narrow strip at 40°N. The 
geostrophic component is of the same sign as the Ekman component, and therefore 
reinforces the SST anomalies set by the Ekman transport. This has been reported in many 
previous studies (Seager et al. 2001; Schneider et al. 2002; Kwon and Deser, 2007) and it 
undermines the earlier notion that the geostrophic adjustment at the KOE is the point of 
origin of the phase reversal of the PDO (Latif and Barnett 1994). The role of Rossby 
waves and delayed response at the KOE were reported to be crucial components of the 
PDO in some model and observational studies (Deser et al. 1999; Seager et al. 2001; 
Kwon and Deser 2007), whereas we find that the SST is largely controlled by the surface 
layer response. It is worth noting that in this study the Ekman transport is sizable even 
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when 10-year low pass filtered wind stress data is used, suggesting that low frequency 
variability in the Aleutian low is the key-controlling factor. 

The lack of preferred periodicity implies that a coupled atmosphere-ocean oscillation is 
unlikely an important mechanism for the PDO simulated in GEOS-5 as opposed to 
CCSM2.0 where a statistically significant 16-year peak was found (Kwon and Deser 
2007). This is further substantiated by the weak atmospheric feedback to the SST (only 
25% of the initial Aleutian low forcing). Additionally, it appears that the atmospheric 
response never matures into the monopole structure of the Aleutian low in order to set up 
a self-sustained oscillation, instead projecting onto the second EOF of the SLP, the NPO. 
The NPO in turn forces an SST pattern, the NPGO, which is resolved as the second EOF 
of the SST. 

In agreement with the previous studies (Chhak et al. 2009), the NPO is an atmospheric 
forcing pattern for the NPGO. However, the PDO and the NPGO, although resolved as 
orthogonal EOFs, are correlated at ~0.4 when the PDO leads the NPGO by 3 years. This 
finding, along with the observation that the PDO’s atmospheric response projects onto 
the NPO, suggests a loose atmospheric link between the PDO and the NPGO. This 
involves an initial Aleutian low anomaly and its oceanic response within 1 year as the 
PDO, and the PDO’s atmospheric response, which matures in 2-3 years, and its reflection 
in the SST in the form of the NPGO. In this scenario, at least in some cases, one could 
expect an NPGO pattern developing after 2-3 years of a PDO peak phase. 
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Previous model studies and results from model inter-comparisons suggest that a 
somewhat realistic spatial structure of the PDO in the North Pacific is a feature 
consistently simulated across the models (Kwon and Deser 2007; Furtado et al. 2011; 
Deser et al 2012). This study provides an additional verification for that and complements 
previous studies. A statistically significant correlation between the PDO and the Aleutian 
low that underlines the atmospheric control over the North Pacific SST appears to be 
another feature consistent among the models. Furtado et al. (2011) notes that this relation 
is captured in the majority of the IPCC AR4 models. Many models, however, fail to 
capture the tropical ENSO connection, overestimate the North Pacific SST variance, and 
differ among the timescale selection. The fact that the PDO is almost entirely of 
midlatitude origin in many models indicates a possible overestimation of midlatitude 
atmospheric control over the decadal SST variability. The enhanced variance over the 
northwestern Pacific would need further attention and may be related to mean biases such 
as too shallow simulated mixed layer as noted by Thomson and Kwon (2010) in 
CCSM3. The overestimation of the low frequency SST variance can also lead to too 
optimistic model predictability estimates. The dominance of Ekman response over the 
geostrophic adjustment indicates predictability of phase reversal of the PDO may be poor 
in this model although persistence of a particular phase may provide some skill. 
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Figure captions 

Figure 1. The (a) DJF seasonal climatological SST (°C) based on 350 years of model 
simulation and (b) difference between model and Hadley Centre SST observation (1870- 
2011). The variance of DJF SST anomalies (°C 2 ) for (c) model simulation (350 years) 
and (d) Hadley Centre SST (1870-2011). 

Figure 2. Spatial pattern of first EOFs 1 and 2 (a, b) in the model and EOFs 1 and 3 (c, d) 
in observations based on PC-SST regression (°C per standard deviation of the PC). The 
EOF is computed from de-trended unfiltered DJF seasonal SST anomalies over the region 
120°E-100°W, 20°N-60°N. The corresponding PCs (grey bars) and their 5-year running 
means (e, f) for the model and (g, h) and observations in standard deviation units. 
Observation is from the Hadley Centre SST (1870-2010). Percentage of explained 
variance (noted in plot labels) is based on total variance values of 0.36 for the model and 
0.21 for the observation. 

Figure 3. Power spectrum of the (a) PCI and (b) PC2 of the SST EOFs (thick black line). 
An average spectrum of 2 141 -year long segments from the 350-year long simulation is 
plotted. Thin solid line is spectrum of best-fit AR1 process and dotted line is associated 
10% confidence level. 

Figure 4. Power spectrum of the (a) PCI and (b) PC3 of EOF modes of observational 
SST for the period 1817-2011. Thin solid line is spectrum of best-fit AR1 process and 
dotted line is associated 10% confidence level. 

Figure 5. Cold and warm composites of monthly SST (shading) and SLP (contour) 
anomalies based on DJF Nino3 index based on (columns 1 and 3, from left) model 
simulation (350 years) and (columns 2 and 4 from left) for Hadley Centre SST and 
NOAA 20 th Century Reanalysis SLP (1871-2010). Months are noted in panel labels and 
the numbers in bracket denote number of months leading or lagging from January. 
Contour levels are +(-) 4, 3, 2, 1.5, 1, 0.5 and 0. Dotted lines represent negative contours. 

Figure 6. Scatter plot of the PCI of the DJF SST anomalies (PDO index) against (a) 
North Pacific Index (NPI) and (b) against the revised NPI that is obtained after removing 
the ENSO signal by linear regression (NPI*). The right panels show the PDO index 
(black curve in c, d) and their respective AR1 forecasts (red curve in c, d) (c) using the 
total NPI and (d) using the NPI*. Correlation coefficients between PDO and the NPI or 
NPI* form the /? parameter for the AR1 forecasts in (c, d). All data are from the model 
simulation. 

Figure 7. EOFs 1 and 2 of DJF (left) SLP and (right) wind stress curl anomalies for the 
model simulation. Percentage of variance (noted in plot labels) is based mean variance of 
8.8 for SLP and 3.9E-5 for wind stress curl. The EOFs are multiplied by the standard 
deviation of the PCs to have units (a, b) hPa and (c,d) Nm" 3 10' 7 . 
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Figure 8. The (top) lead/lag correlation between the (open circle) PCI and PC2 of SST 
and (closed circle) PCI and PC2 of wind stress curl. The (bottom) Lead/lag correlation 
between (open circle) PCI of SST and PCI of wind stress curl, (closed circle) PC2 of 
SST and PC2 of wind stress curl and (open square) PCI of SST and PC2 of wind stress 
curl. Negative lags in the top panel indicate PCI leads PC2 and in the bottom panel 
atmosphere leads SST. Horizontal lines with a symbol represent 5% significance level for 
based on two-sided t-test. Since only positive correlation coefficients are discussed, 
significance level is noted only for that. All time series are 6-year low pass filtered. All 
data are from the model simulation. 

Figure 9. Lead/Lag regressions (shading) and correlations (contours) of (column 1) SST 
(°C), (column 2) SLP (hPa) and (column 3) wind stress curl (10 s Nm" 3 ) w.r.t the 
normalized PCI of SST (i.e., PDO index). Lags are in years. Negative (positive) lag 
indicates PDO lags (leads) the field. All variables are from the model simulation and are 
6-year low pass filtered. Hatching indicates significance at 5% level from a two-sided t- 
test. 

Figure 10. Same as Fig. 9: but w.r.t the PC2 of the SST or the NPGO index. 

Figure 11. Lead/lag regressions (shading) and correlations (contours) of the net surface 
heat flux (Wm 2 per °C) with respect to the KOE index. Positive net heat flux is directed 
into the ocean. Hatching indicates 5% significance level based on a two-sided t-test. 

Figure 12. Simultaneous regressions (shading) and correlations (contours) of temperature 
advection by (left column) Ekman and (right column) geostrophic currents in the mixed 
layer (Wm 2 per °C) with respect to the KOE index. The anomalous transport is 
decomposed into (a, e) advection of temperature anomaly by mean currents or V ■ VT', 
(b, f) advection of mean temperature by anomalous currents or V' ■ VT, (c, g) advection of 
temperature anomaly by anomalous currents or V' ■ VT', and (d, h) sum of all three or 
total anomalous transport, where V stands for total Ekman or geostrophic velocity vector 
in the mixed layer. The overbar denotes climatological mean and prime denotes deviation 
from climatology. Hatching indicates 5% significance level. 
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Figure 2. Spatial pattern of first EOFs 1 and 2 (a, b) in the model and EOFs 1 and 3 (c, d) 
in observations based on PC-SST regression (°C per standard deviation of the PC). The 
EOF is computed from de-trended unfiltered DJF seasonal SST anomalies over the region 
120°E-100°W, 20°N-60°N. The corresponding PCs (grey bars) and their 5-year running 
means (e, f) for the model and (g, h) and observations in standard deviation units. 
Observation is from the Hadley Centre SST (1870-2010). Percentage of explained 
variance (noted in plot labels) is based on total variance values of 0.36 for the model and 
0.21 for the observation. 
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970 Figure 3. Power spectrum of the (a) PCI and (b) PC2 of the SST EOFs (thick black line). 

971 An average spectrum of 2 141 -year long segments from the 350-year long simulation is 

972 plotted. Thin solid line is spectrum of best-fit AR1 process and dotted line is associated 

973 10% confidence level. 
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977 Figure 4. Power spectrum of the (a) PCI and (b) PC3 of EOF modes of observational 

978 SST for the period 1817-2011. Thin solid line is spectrum of best- fit AR1 process and 

979 dotted line is associated 10% confidence level. 
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Figure 5. Cold and warm composites of monthly SST (shading) and SLP (contour) 
anomalies based on DJF Nino3 index based on (columns 1 and 3, from left) model 
simulation (350 years) and (columns 2 and 4 from left) for Hadley Centre SST and 
NOAA 20 th Century Reanalysis SLP (1871-2010). Months are noted in panel labels and 
the numbers in bracket denote number of months leading or lagging from January. 
Contour levels are +(-) 4, 3, 2, 1.5, 1, 0.5 and 0. Dotted lines represent negative contours. 
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Figure 6. Scatter plot of the PCI of the DJF SST anomalies (PDO index) against (a) 
North Pacific Index (NPI) and (b) against the revised NPI that is obtained after removing 
the ENSO signal by linear regression (NPI*). The right panels show the PDO index 
(black curve in c, d) and their respective AR1 forecasts (red curve in c, d) (c) using the 
total NPI and (d) using the NPI*. Correlation coefficients between PDO and the NPI or 
NPI* form the /? parameter for the AR1 forecasts in (c, d). All data are from the model 
simulation. 
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1002 Figure 7. EOFs 1 and 2 of DJF (left) SLP and (right) wind stress curl anomalies for the 

1003 model simulation. Percentage of variance (noted in plot labels) is based mean variance of 

1004 8.8 for SLP and 3.9E-5 for wind stress curl. The EOFs are multiplied by the standard 

1005 deviation of the PCs to have units (a, b) hPa and (c,d) Nm' 3 10' 7 . 
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Figure 8. The (top) lead/lag correlation between the (open circle) PCI and PC2 of SST 
and (closed circle) PC 1 and PC2 of wind stress curl. The (bottom) Lead/lag correlation 
between (open circle) PCI of SST and PCI of wind stress curl, (closed circle) PC2 of 
SST and PC2 of wind stress curl and (open square) PCI of SST and PC2 of wind stress 
curl. Negative lags in the top panel indicate PCI leads PC2 and in the bottom panel 
atmosphere leads SST. Horizontal lines with a symbol represent 5% significance level for 
based on two-sided t-test. Since only positive correlation coefficients are discussed, 
significance level is noted only for that. All time series are 6-year low pass filtered. All 
data are from the model simulation. 
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1022 Figure 9. Lead/Lag regressions (shading) and correlations (contours) of (column 1) SST 

1023 (°C), (column 2) SLP (hPa) and (column 3) wind stress curl (10‘ 8 Nnf 3 ) w.r.t the 

1024 normalized PCI of SST (i.e., PDO index). Lags are in years. Negative (positive) lag 

1025 indicates PDO lags (leads) the field. All variables are from the model simulation and are 

1026 6-year low pass filtered. Hatching indicates significance at 5% level from a two-sided t- 

1027 test. 
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1029 

1030 Figure 10. Same as Fig. 9: but w.r.t the PC2 of the SST or the NPGO index. 
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1034 Figure 11. Lead/lag regressions (shading) and correlations (contours) of the net surface 

1035 heat flux (Wm’ 2 per °C) with respect to the KOE index. Positive net heat flux is directed 

1036 into the ocean. Hatching indicates 5% significance level based on a two-sided t-test. 
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Figure 12. Simultaneous regressions (shading) and correlations (contours) of temperature 
advection by (left column) Ekman and (right column) geostrophic currents in the mixed 
layer (Win 2 per °C) with respect to the KOE index. The anomalous transport is 
decomposed into (a, e) advection of temperature anomaly by mean currents or V ■ VT', 
(b, f) advection of mean temperature by anomalous currents or V' ■ VT, (c, g) advection of 
temperature anomaly by anomalous currents or V' ■ VT', and (d, h) sum of all three or 
total anomalous transport, where V stands for total Ekman or geostrophic velocity vector 
in the mixed layer. The overbar denotes climatological mean and prime denotes deviation 
from climatology. Hatching indicates 5% significance level. 


46 




